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Abstract 

The renormalized mass of the bottom quark is calculated at the two loop level 
to order (9(a s GVM t 2 ) in the MS renormalization scheme. Different strategies 
for the computation are outlined. The result is applied to the partial decay 
rate T(H — > bb) of the Higgs boson into bottom quarks. Expressing the width 
in terms of the running mass instead of the bottom pole mass allows to treat 
the 0(a s GFMf) radiative corrections on the same footing as is commonly 
used in pure QCD calculations. The numerical values for the corrections are 
given and the sizes of different contributions are compared. 
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1 Introduction 



Studying the properties of the Higgs boson, once it is discovered in future particle ac- 
celerators, will be the prime tool to experimentally probe the details of the electroweak 
symmetry breaking mechanism in the Standard Model. Of particular interest will be the 
Higgs boson decay into bottom quarks, since the decay channel H — > bb dominates in the 
intermediate Higgs mass range Mh < 2M^. This process will be even more important, 
if possible hints for new physics effects in the reported descrepancy [|TJ between the mea- 
sured partial Z boson width Rf, into bottom quarks and its theoretical prediction should 
happen to substantiate. Similar effects might then also be visible in Higgs decays H —>■ bb 
and emphasize the need for precise Standard Model predictions to T(H — > bb). 

As a consequence much work has been spent on the calculation of radiative correc- 
tions to Higgs processes in the past and excellent reviews on Higgs phenomenology can 
be found in the literature |2|, [|. Previous works concerning the partial rate T(H — > bb) 
include electroweak one loop corrections (|, ||, |J , the calculations of universal and nonuni- 
versal corrections of the order 0{a s GpM^) [0, |8|, [9|, PH , and recently even a three loop 
0(a 2 s GFMt) calculation was presented |TT|, [T^j. Nonuniversal corrections to the vertex 
Hbb involve the virtual top quark through Higgs ghost exchange. Their top mass en- 
hancement oc rn? t due to Yukawa couplings distinguishes them from similar vertices of the 
Higgs boson to other quark flavours. 

In our earlier work || the diagrams of Figure [l] were considered in the heavy top limit 
M 4 2 3> Mfj. The two loop (!?(a s GpM t 2 ) relation between the bare mass m and the on- 
shell (OS) mass Mb of the bottom quark was presented and the corrections to the partial 
Higgs width were expressed in terms of M^. 

The MS renormalization scheme on the other hand is the commonly used renormaliza- 
tion prescription in higher order QCD calculations. Apart from calculational convenience 
its concept of the running bottom mass allows the absorption of large logarithms 
ln(M^ /Mfj) (see e.g. [TJ], [TJ3|, |TB[) and causes the perturbation series to converge more 
rapidly than in the OS scheme. It is therefore of obvious interest to adopt the notion 
of the running MS mass mj in the Higgs decay rate also for the case when electroweak 
corrections are included. 

For this reason we have calculated the two loop relation of order 0{a s GpM 2 ) between 
the on-shell mass and the MS renormalized mass of the bottom quark (for a discussion at 
the one loop level see [17]]). This transformation from one renormalization scheme to the 
other allows to express T(H — > bb) to the order 0(a s GfM, 2 ) throughout in terms of the 
running mass rhf,. 

The problem is approached in four different ways. All methods are leading to the same 
answer and thus provide powerful crosschecks beyond the standard consistency checks such 
as gauge invariance. 

In order to introduce our notation let us start from the bare Lagrangian and consider 
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the bare fermion propagator for the bottom quark 

S, 1 = i{m -i>- m S° 5 - (1) 

We have not written the term ^7sE^ for notational simplicity, since for all quark mass 
relations below becomes relevant only in higher order electroweak corrections 0(G%), 
which we do not consider in this work. We quote our previous result for the 0{a s GpMf) 
bottom pole mass 

M b = ™o^ff (2) 

in Eq.(^6|) of the appendix, where for later convenience and are expressed in terms 
of the MS masses ffi b and mj. 

By rescaling its parameters the bare Lagrangian can be written as the sum of the 
renormalized Lagrangian and the counterterm Lagrangian. Our interest focuses on the 
renormalization constants Z 2 and Z m relating the bare wavefunction and mass of the 
bottom quark to their renormalized equivalents 

= Zl /2 ^, m = Z m m b . (3) 

Here we adopt the MS renormalization scheme as is indicated through bars. The renor- 
malized bottom quark propagator accordingly reads 

Sft 1 = Z 2 S 1 

(4) 

= i(m b ~i>- rhbTjs - f&v + (Z 2 Z m - l)rh b - i>{Z 2 - 1 

For the determination of the MS bottom mass we perform our calculations according to 
the following different strategies. 

In Section 2.1 the overall counterterms Eg" 7 ", £y- T to the bottom selfenergy are com- 
puted in the MS scheme. With 



1-Zy 



CT 

(5) 



Z 2 Z m — 1 + Ti S T 

one obtains the relation between the MS and bare masses 

1 - fP T 

m b = mo i + fCT - ( 6 ) 

In combination with Eq.(|2|) this leads to the transformation rule between OS- and MS 
masses of the bottom quark. 



M - a + Sg T )(l-£°; m 

Mb = mfc (i-s^)(i-s°v ( 7 ) 
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In Section 2.2 a different approach is used to verify the findings of Section 2.1. The 
renormalized bottom quark propagator Eq.(|]) is rewritten in the form 

Sj, 1 = i [m b (l - £ s + Sg T ) - i (l + Z v - S£ T ) } (8) 

with Eg = Z2Z m TP Sl XV = ZzEfy. We check by explicit calculation of the finite parts of 
the bottom quarks self energies Eg n , Ey n that the relation 

[I - £a + S? T ) 



M 6 = m 6 - 



m 6 



E 



V I 



;i+s{n 



(9) 



is indeed equivalent to the prescription Eq.(|7j). 

In Section 2.3 our problem is considered from a third point of view, which becomes 
transparent by expressing the renormalized quark propagator in the following form 

Sk 1 = iZ 2 (Z m m b -jf - Z m rh b Y? s -^S^). (10) 

This expression is finite if the bare parameters in TP SV are substituted in favour of the 
renormalized ones. One therefore can solve for Z 2 and Z m recursively, i.e. loop by loop. 
The renormalization constant Z m leads then to the same result for fh b as in the previous 
sections. 

Finally we demonstrate in Section 2.4 that another simple derivation of the result is 
possible, based on the earlier determination of the bottom pole mass and leading to the 
same fh b again. 

The results are then applied in Section 3 to the partial decay rate T(H — ► bb). The 
numerical size of the corrections are given and the renormalization scheme dependence is 
discussed. 



2 Calculation of the MS Renormalized Bottom Mass 
2.1 Approach 1: Counterterms 

We calculate the MS counterterms on a graph by graph basis in this section. Per definition 
of the MS scheme counterterm vertices consist of pole terms only and are therefore easier 
to compute than full diagrams. The integrals represented by the graphs in Figure [I] involve 
several mass scales. Via electroweak interactions the top quark and the Higgs ghost come 
into play with their respective scales M t and Mw- In the heavy top limit — » oo one 
has M t 2 . Since we consider only the leading term oc M f 2 in the power series of the 

inverse top mass, one can neglect Mw right from the beginning. As a consequence the 



electroweak gauge parameter drops out trivially. The heavy mass expansion |L8], [19], gy, £L 
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has developed into a well established technique and was sucessfully used in a number of 



applications. For a more detailed description the reader is referred for example to f22 |. 
The main virtue of this method is the factorization of a multiloop integral containing the 
heavy top quark into an integral with less number of loops and massive tadpole integrals. 

This decomposition is operative in our problem as well. Two loop integrals eventually 
factorize into a one loop tadpole and a one loop propagator integral, where the latter 
involves two scales, namely the bottom mass and the external momentum. However, being 
interested in the pole parts only, the matter simplifies even more. Since the pole terms 
are independent of masses and momenta, one can conveniently nullify either of them. 
Care must be taken that no spurious infrared divergencies are introduced in this way. In 
our case we have obtained the pole parts to £5 by setting the external momentum to 
zero, thus reducing the massive propagator integral to a tadpole integral. Similarly, for 
the computation of £y the bottom mass is nullified. The resulting massless propagator 
integral is conveniently computed with the help of MINCER |S| which is based on the 
symbolic manipulation program FORM |[24|| . 

The counterterms of the one loop diagrams "QCD" and "EW" of Figure [l] are simply 
given by their pole terms obtained in the above described manner. On the two loop level 
the situation is somewhat more involved, since the diagrams "IN" , "OUT" and "LEFT" 
contain ultraviolet divergent subgraphs. As is indicated in Figure ||, these sub divergences 
have to be subtracted in order to arrive at the overall divergence of the corresponding 
diagrams. The removal of the sub divergences results in local counterterm vertices, which 
we list in the appendix. It can be seen that indeed all logarithms have dropped out. 

Whereas the counterterms are still gauge dependent, the QCD gauge parameter £ s 
cancels in the following expression for the bottom mass: 



m b = m 



1 + ^s T (11) 

f, , at s 1 _ 3 1 a s _ 2 

= TO Q H Vx t --A x t - 

V 7T e 2 e 7r e 

with x t = G Frhf / S\/2tt 2 . This leads to the transformation between the pole and the MS 
mass of the bottom quark 



M b = m b 



(1 + Eg}(l-E|) 

:i-^ T )(i-4) 

'5 3 u 2 




m b {1 -\ ( - + In — } + x t \ - + - In _ 2 



4 2 fhfj 

a s f 9 ,, . 5 , u 2 5 , u 2 3 , , u 2 3 , u 2 , u 2 
^42- 4C(2) + 2 ln ^n ln £ 2 + 2 ln ^ 2 + 2 ln ^ ln | 2 



(12) 
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2.2 Approach 2: Finite parts 



As a cross check of the result Eq.(|12D we now want to recalculate it in a different way, 
namely by employing only the finite parts of the corresponding bottom self energy graph 
as given in Eq.((^). The finite part of a diagram 

yifin yifull spsub yiCT /iq\ 

n s,v — ^s,v ^s,v ^sy 

is obtained by subtracting the overall counterterm S^y and the counterterm with the 
subdivergence Egy from the full diagram figy. Pictorially this procedure is visualized in 
Figure [5]. Notice that £gy contains both pole and finite terms. One therefore cannot use 
the nullification procedure of the previous section to simplify the calculation. 

Instead it is possible to simplify integrals by evaluating them on the mass shell p 2 = Mf 
25 1 using the expansion 



b 



^v(f)= Sg?(l) + (f ~ l) EWl) 



/ (14) 
Sg?(l) + 2(S S + Sv) T? 8>V {1). 



The derivatives £ sy = dT^sy /d{fh 2 /p 2 ) may be conveniently obtained through deriva- 
tions with respect to rhb, thus raising the power in the denominator of the integrand. This 
procedure may also be applied for the calculation of subdivergence counterterms, where 
the corresponding expansion reads 



ml 



ftfv [f) = + 2 + t™) t' sy {l). (15) 

The expressions for the finite parts of the various contributions are listed in the ap- 
pendix. They lead to the relation between pole and MS bottom mass 



M b = m b 



;i - sg n ) 



2 




a s ( 4 1 _ 

b 




5 /i 2 5 /i 2 3 2 /i 2 3 /i 2 /i 2 

We find agreement with Eq.(|12|) 



4C(2) + ^ln^ + ^ln^ + ^ln 2 ^ + ^ln^ln _ , 



(16) 



2.3 Approach 3: Renormalization Constants 

In our third method we proceed along a path which deals directly with the renormalization 
constants Z 2 and Z m . To explain how both Z 2l Z m are computed iteratively loop by loop, 
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it is convenient to consider the the renormalized fermion propagator in the following form: 

Sx 1 = iZ 2 (z m m b - jl - Z m rh b Y? s - #fE y ) (17) 



Here the bare bottom selfenergies TP SV = S^y (m^o, mt,o) + ^sv ( m &,o, m t,o) receive con- 
tributions from the one and two loop diagrams of Figure |l|. The explixcit arguments shall 
emphasize that all parameters are understood as bare quantities. In general the parame- 
terlist would also include coupling constants, gauge parameters etc. If we now substitute 
the bare masses in favour of their renormalized counterparts at a given loop level, the 
functional form of the selfenergies does not change in that given order, but additional 
contributions of higher order are induced: 

4!y°Koi m *,o) = £${?(m 6 , ™*) + S S/ind ( 18 ) 

Let us first consider the one loop case. Having expressed Eq.(|l7]) entirely in terms of 
renormalized quantities, the renormalization constants Z 2 ,Z m must be such that the in- 
verse quark propagator is finite and, stated more precisely for the MS scheme, that the 
poles cancel. According to the Lorentz structure this results in two equations 



m b Z 2 Z m ( 1 - S^ 1)0 (m fc , m t ) ) = finite 

(19) 

l + 4 1)0 (m 6 ,m t )) = finite, 
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which can be solved for Z 2 and Z m . The solution for Z m leads to the one loop result for 
the MS mass fh b . 

The procedure can then be repeated for the two loop case. Besides the two loop result 
£^ y also the induced second order terms E^ 2 yind from the transition to the renormalized 
parameters at the one loop iteration have to be taken into account. Solving the corre- 
sponding system of equations gives the renormalization constants at the two loop level. 
Inversion of Eq.(|3|) with 

Z m = 1 - x t - 1 x t — . (20) 

2e 7r e 7r \e z ej 

indeed confirms the result Eq. ([□])• 

2.4 Approach 4: Derivation from OS Mass 

Having approached the problem from three different sides, let us demonstrate, how the 
MS bottom mass can be derived in another elegant manner. We start with the following 
ansatz for the relation between the bare mass and the MS mass of the bottom quark (a 
similar method was used in |26||) and insert it into Eq.(0): 



M b = m b \ 1 + — - + x ttQ - + —x tfi - + - )\ — — £ 21 
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The bare top mass in x t)0 is substituted through the renormalized MS mass and the terms 
following the curly bracket are taken from the pole mass calculation in Eq. ( |26| ) . 

The crucial step is to require that the pole mass Mb on the LHS as a physical quantity 
must be finite. This translates into the requirement that all coefficients of 1/e poles on the 
RHS must vanish. Thus one obtains four equations which can be solved for the unknown 
coefficients a, b, c, d. Two additional equalities follow from the fact that the logarithms of 
the pole terms cancel separately and serve as a consistency check for the solutions 

a = -l, b=-^, c = 0, d = -2. (22) 

Insertion into Eq.(pl|) produces again the result Eq.(|l2"l) and Eq.flTBf). 



3 Application to the Higgs Decay H —> bb 

In this section we apply our result to the partial Higgs boson decay rate || |10 

Y{H - bb) = T M 2 jl + X t + ^X t (-1 - 4C(2) - 2 In j (23) 

where T = 3GfM h /4\/2ii, X t = GpM 2 /S\f2-K 2 and the renormalization scale is chosen 
as /i 2 = Mfj. For our following numerical discussion we use as input values a bottom 

(5) 

pole mass of M& = 4.7 GeV and a top mass of M t = 176 GeV. Based on Aq CD = 233 
MeV the running strong coupling constant ranges between a s (M H = 70 GeV) = 0.125 
and a s (Mn = 130 GeV) = 0.114 where a s (/i) is defined for five active flavours. We now 
express the above formula for the width in terms of MS masses and obtain 



^ J 1 , - ,7 . ..Ml' 



T(H -> bb) = T mi <1 + x t - + 3 In 



2 m 




M 2 

12t(2) + 26 In— f + 3 In 

mi 




(24) 



Notice that the transformation of Eq.(|23|) implies that the first order QCD corrections 
T M 2 (a s /7r)[3 — 2ln(M^/M b 2 )] give rise to a contribution of order 0(a s G F rh 2 ) as well. 
Based on the given values for the pole masses the corresponding running masses amount 
to m b = 2.84/2.69 GeV and m t = 179.44/170.04 GeV for M H = 70/130 GeV. In Figure | 
we plot the corrections of orders Gpm\ and Gpm\ + a s Gpm 2 according to Eqs. (|23|) and 
(|4]). The curves are strongly characterized by the linear rise in M# due to the overall 
factor. For the on-shell result the QCD screening of the leading electroweak corrections is 
clearly visible. The MS curves indicate that the two loop contribution is less important 
than for the OS scheme and suggest a better convergence behaviour of the perturbation 
scries. 
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An inspection of Eq.(|24D reveals that all large logarithms ln(Mfi/M^) have dropped 
out. Their absorption into the running bottom mass favours the use of the MS mass 
over the OS-mass. There is no such strong preferance with respect to the top mass, 
considering that the scales of the Higgs and the top are not as far apart as the Higgs and 
the bottom scales. Corresponding logarithms ln(Mj^/M t 2 ) therefore cannot be considered 
as particularly dangerous. Instead one might tend to use the top pole mass as a quantity 
which is by definition more feasible in experiments. In this case the Higgs decay rate can 
be rewritten into the following form 

T(H^bb) = T ml (l + ^+X t (Z + 3m^f| 




Mi. , 

12<(2 ) + 17 ln^-3^f' 

30.717 -HAlm* M ^ 



(25) 



3 M} 9 m\ 



Several groups have contributed to the calculation of QCD corrections which we have 
included in the formula in first JT3], [P], ^7], ^| and massless second order |5C 



Quadratic bottom mass corrections in second order [161, HI anc ^ top quark contributions 



|32|, [33], j34j are also available. We have not written these pieces into Eq.(^), but included 
them in our numerical analysis. 

As was noticed in |33| the logarithms c? s ln 2 (M^-/m^) originate from flavour singlet type 



diagrams. They are not present in the rate for the decay into hadrons, but are introduced 
when the pure gluonic channel is subtracted. In Figure |5| the contributions coming from 
the orders a s , a 2 ^ Gi7M t 2 , a s GVM t 2 are compared, normalized to the Born term TBom = 
r m 2 . The electroweak corrections may carry different sign depending on the Higgs mass. 
However, as compared to the QCD corrections, where the first order contributes about 
20% and the second order about 5% to the corrections, the electroweak contributions are 
small. With 0(G F M?) = -6.6/5.4 per mille and 0(a s G F M?) = -4.3/ - 3.9 per mille 
for Mu — 70/130 GeV these effects become relevant for high precision experiments at the 
percent level. 
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A Appendix 



The result from || for the relation between the pole mass and the bare mass of the bottom 
quark according to Eq.(§) is reproduced below. For convenient use in Section 2 the bare 
masses in £5 and are already transformed into the MS renormalized ones. 



M h 



mo- 



a s 1 1 4 , a 2 
1 + — - + -+ln-E^ + e 



7r \ e 3 ml 



^ + ic(2) + ^ln^ + iln 2 



/r 



mi 



■v-,2 



,31 5 3 , u 2 
~^~ x t [ 7: — h — + — In ——y + e 



2e 



H Zt 

7T 



^1 

+ 



4 

21 



mi 

3, /i 2 
h - In At 

4 2 m 2 t 



- + -C(2) 
8 4 SV ; 



5, l? 

4 



3, /i 2 
2 



2 4 m£ 



9, 2 /i 
+ - In 2 - 

4 



m 2 



+ -ln 2/i 



13, /i 2 

- 2 + T ln ^2 
4 m 2 

+ ^ln4ln^ 

it 



4 m 2 



(26) 



The overall counterterms are given for the different diagrams: 



(QCD) 



Y,y T (QCD) 



7r e V 
a s 11 
7r e 3 



(27) 



E£ T (£W0 



2e 



(IN) 



S^ T (IN) 



a* 



7T 

a s _ 
— £t 

7T 



e 2 12^ s + e 24 ^ £ 



3 3*' 



(28) 



(OUT) 



££ T (OUT) 



• <4 (~ + ^ 



a 

7T 

a: 



1 

6' 



7T 



1_J_ 1 

F 12^ + 7 



6 6^ 



-06 



(29) 



Eg T (LEFT) 
Sy T (LEFT) 



/in 2 3 e 



1 

e 
1 

+ - 

e 



1 1 

■2 + 5«- 



1 J_ 

3 + 12^ 



(30) 
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The finite parts of the one and two loop contributions read as follows: 



f.r, { -2 -- 2 In 



S^(l-loop)= ^{^ s + kln^ 
7r 3 3 ?7r 



■ 3 1, /i 

4 2 mf 



6 

2 ~1 



6 3"V""m? 



1 i 2 ^ , 1 > i V 2 i 
-In — + -£ s ln — In — 

2 6 rrif ml 
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Figure Captions 

Figure 1: Order 0{a s GpM^) self energy diagrams for the bottom quark. Thin line: 
bottom quark, thick line: top quark, curly line: gluon, dashed line: Higgs ghost. 

Figure 2: Counterterm diagrams up to order 0(q s GfM ( 2 ), 

Figure 3: Finite terms for diagrams up to order 0{a s GpMf). 

Figure 4: Corrections to T(H — > bb) in terms of pole masses (upper curves) and MS 
masses (lower curves). The solid lines are the O(GFmf) contributions and the dashed 
lines the sum of O(GFm^) and 0(a s GFmf). 

Figure 5: Corrections T(H — > bb) separately for the orders 0(a s ) (solid curve), 0(a 2 s ) 
(dashed-dotted curve), 0(GVM t 2 ) (dashed curve) and 0(a s G F M?) (dotted curve . 
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